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ABSTRACT 

We examine a sample of 45 Mg II absorption-selected systems over the redshift range 0.4 to 1 .4 in order to 
better understand the range of physical conditions present in the interstellar and halo gas associated with inter- 
mediate redshift galaxies. Mgll and Fell absorption profiles were observed at a resolution of ~ 6 km s" 1 with 
HIRES/Keck. Lyo and Civ data were measured in FOS spectra obtained from the Hubble Space Telescope 
archive (resolution ~ 230 km s _1 ). We perform a multivariate analysis of W r (Mgll), ^(Fell), W r (Clv) and 
W r (Ly a) (rest-frame equivalent widths) and the Mg II kinematic spread. There is a large range of high-to-low 
ionization properties and kinematics in intermediate redshift absorbers, that we find can be organized into five 
categories: "Classic", "C iv-deficient", "Single/Weak", "Double", and "DLA/Hl-Rich". These categories arise, 
in part, because there is a strong connection between low-ionization kinematics and the location of an absorber 
on the W r (Clv)-W r (Mgll) plane. Using photoionization modeling, we infer that in most absorbers a significant 
fraction of the Civ arises in a phase separate from that giving rise to the Mgll. We show that many of the Civ 
profiles are resolved in the FOS spectra due to velocity structure in the Civ gas. For 16 systems, the galaxy 
Mk, Mb, B — K, and impact parameters are measured. We compare the available absorption-line properties 
(taken from Churchill et al. 1999, Paper I) to the galaxy properties but find no significant (greater than 3 a) cor- 
relations, although several suggestive trends are apparent. We compare the locations of our intermediate redshift 
absorbers on the W r (C lv)—W r (Mg n) plane with those of lower and higher redshift data taken from the literature 
and find evidence for evolution that is connected with the Mg II kinematics seen in HIRES/Keck profiles of Mg II 
at z > 1.4. We discuss the potential of using the above categorizations of absorbers to understand the evolution 
in the underlying physical processes giving rise to the gas and governing its ionization phases and kinematics. We 
also discuss how the observed absorbing gas evolution has interesting parallels with scenarios of galaxy evolution 
in which mergers and the accretion of "proto-galactic clumps" govern the gas physics and provide reservoirs for 
elevated star formation rates at high redshift. At intermediate and lower redshifts, the galaxy gaseous components 
and star formation rates may become interdependent and self-regulatory such that, at z < 1, the kinematics and 
balance of high and low ionization gas may be related to the presence of star forming regions in the host galaxy. 

Subject headings: quasars — absorption lines; galaxies — evolution; galaxies — halos 



1. INTRODUCTION 

It is well known that individual galaxies fall into a fairly 
clear-cut morphological classification scheme, as first proposed 
by Hubble (1936). Over the decades, Hubble's scheme has been 
expanded upon and exploited for quantifying galactic evolution, 
galactic stellar populations and star formation efficiencies, the 



relative gas-phase component to stellar component in galaxies, 
and the effects of environment on galaxy formation. Similarly, 
one of the central motivations for studying intervening quasar 
absorption lines, especially those selected by the presence of 
metal lines, is that they also provide insights into galactic evo- 
lution, not of the stars and stellar dynamics, but of the chemi- 
cal, ionization, and kinematic conditions of interstellar and halo 



1 Based in part on observations obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among Caltech, the University of California, 
and NASA. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. 

2 Based in part on observations obtained with the NASA/ESA Hubble Space Telescope, which is operated by the STScI for the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS5-26555. 

3 Visiting Astronomer at the W. M. Keck Observatory 

4 Center for Gravitational Physics and Geometry 

5 NSF Young Investigator 



2 



Mgll ABSORPTION SELECTED GALAXIES 



gas. 

Historically, absorption line systems have been classified 
a taxonomic system by their gas cross sections. Categ< 
by increasing Hi column densities are the Lya forest cl 
with AT (Hi) < 10 15 cm -2 , the sub-Lyman limit systems 
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~ 2 , the Lyman limit break systems 
2 , and the damped Lya systems, 
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A (Hi) > lCT " cm 

ing AT (Hi) > 10 20 3 cm~ 2 . Metal-line systems are eithir 
lected by the presence of strong Mg II AA2796, 2803 absoi pt: 
or strong Civ AA1548, 1550 absorption. 

What has not been explored, however, is the taxonomy 
sorption line systems when the H I, Mg II, and C IV absoijpt 
strengths and the gas kinematics are equally considered 
mately, categorizing the relationships between several differ 
absorption properties may provide clues central to understand 
ing the different physical natures of the various types ol 
tems. 

Furthermore, while our empirical knowledge of the 
tion strengths, kinematics, and physical extent of galaxies 
lected by M g II absorption has steadily progressed over th 
decad e (e.g. Bergeron & Boisse 1991 ; panzetta & Bowen" ! 990 



ab- 
ion 
ind- 
ent 



sys- 
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1992; |Bergeron et al. 1992 



^e Brun et al. 19931; Steidel let al 
1994; phurchill et al. 1996| ; puillemin & Bergeron 1997 ), lit 
tie is directly known about Hi absorption and higher ionization 
absorption (esp. Si IV, Civ, Nv, and O VI) in these galaxies 

Arguably, the Mg Il-selected systems are ideally suited 
taxonomic study of absorption systems because: 
(1) those with W r (Mgll) > 0.3 A are known to be 



associated with galaxies ( Bergeron & Boisse 1991 
Dickinson, & Persson 1994; Churchill et al. 1996? and/or 



Yaniy 



galactic metal-enriched environments ([Yanny 1992| 
York 1992). HST imaging has revealed that these galaxies 
a wide variety of line-of-sight orientations and " normal" 
phologies (see |Steidel et al. 1997[ |Steidel 199$ ). Since 
nesium is an a-group element yielded by Type II 
it is expected that the association with galaxies will hold t 
highest redshifts. 

(2) they arise in structures having a five decade 
of Hi column densitie s, including sub-Lyman limit systems 
flChurchill et al. 1999a] ; C hurchill et al. 2000, here after Pat 
Lyman limit systems (e.g. Steidel & Sargent 1992 1 Paper I 
damped Lya systems (e.g! [Le Brun et al. 1997t R ao & 
shek 1998; [Boisse et al. 1998[ ) which means that a large 
of galactic environments will be sampled. 

(3) for z < 2.2 their statistical properties and for z 
their kinematic properties have been th oroughly documented 
(e.g. Lanzetta, Turnshek, & Wolfe 1987[ S argent, Bokserjber 
& St eidel 1988; gteidel & Sargent 19921 Petitj ean & Berj 
1 990; [Steidel & Sargent 1992fc |Churchill 19971; C hurchill 
1999a) which means that the low redshift database is 
place and can be used for evolution studies when higher redshift 
data are obtained. 

(4) they are seen to give rise to a range of ClV and 
high ionization absorption in UV space-based spectra ( 
eron et al. 1994; Paper I), which means that the more gejneral 
ionization con ditions can be studied in detail (e.g. C hurcflill & 
Charlton 1999). 

In this paper, we study the combined available data on a sam- 
ple of 45 Mg II absorption-selected systems at redshifts ~ 0.4- 
1.4. We focus on the Mgll and Fell absorption strength; and 
kinematics, the Lya and Civ absorption strengths (Papjr I), 
and the ground-based derived galaxy luminosities, colors , and 
impact parameters. We perform a multivariate analysis (pabu 
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o objectively quantify systematic differences 
r^upings of the high ionization and neutral hy- 
properties of Mgll selected systems, 
ganized as follows: In § ||, we describe the 
ielection, and the data analysis. In § ||, we 
in the high and low ionization proper- 
e redshift Mg II absorbers using multivariate 
a taxonomy that serves as an objective 
g the variations in Mgll absorber proper- 
estigate the overall ionization and kinematic 
we investigate the relationship between the 
and host galaxy properties and in § || of- 
discussion on the relationships between the 
nd their possible evolution. In § ^ we sum- 
of this work. 



p aints 1 



2. THE DA' A: SAMPLE SELECTION AND ANALYSIS 



We targeted quasar 
ate redshifts that 
Mgn AA2796,28 



absorption line systems at intermedi- 
have been discovered by the presence of a 
)3 doublet (e.g. Lanzetta, T urnshek, & Wolfe 
1987 ; ISargeritTBcksenberg, & Steidel 19881; Steidel & Sargent 



1992). We also in 
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lude in our study the "weak" systems [those 
; 0. 3 A], which are more n umerous in their 
hurchill et al. 1999ah. A detailed ac- 
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these systems, additional 
was available in 
piled by the HST 



e selection is given in Paper I, but we briefly 
properties here. 

rber s were selected from a high resolution 
) survey (phurchill 1997| ) using the HIRES spec- 
al. 1994| ) on the Keck I telescope. For 45 of 
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Johnson & Wichern 1992) on the absorption 



wavelength coverage in the ultraviolet 
HST/FOS archive and the database com- 
QSO Absorption Line Key Project (B ahcall 
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database of Mg II absorbing systems and their 
lent widths are listed in Table [I] The first 
the quasar name, the absorber redshift, and 
frame equivalent width, W^Mgll). The 
columns list the Mg II kinematic spread, ui v (see 
kinematic composition, i.e. the number of 
components (see § |2.2[ ). Columns six, seven, and 
Civ, and Lya rest-frame equivalent widths. 
iv0s the various subsamples used for our anal- 
lists the taxonomic "class" of each system 
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t le data analysis are given in Paper I, includ- 
fitting, the line finding, the equivalent width 
establishment of a redshift zero point for the 
procedure for line identifications, and the tech- 
1 measure Lyman limit breaks, when present. 
Fi jure 2 of Paper I, the Mg II equivalent width 
con iistent with th at of an unbiased sampl e of ab- 
gll) < 1.3 A ( phurchill et al. 1999a[ ). There 
w ith W r (Mg II) > 1.3 A, which may appear as 
stro lg systems. When this may be a concern for 
c iscuss this possibility, 
letection sensitivity, the sample is 72% com- 
II) rest-frame detection threshold of 0.02 A 
to a . 03 A threshold (Figure 3 of Paper I) . In 
Ta ble |], we have designated those systems hav- 
deqection thresholds greater than 0.03 A as Sam- 
sorbers, unresolved Mg II absorption features 
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with W r (MgIl) < 0.03 A cannot be detected. The other sub- 
sample designations, "CA" and IC", are explained in §§ || and 
||, respectively. 

2.1. Galaxy Sample 

For roughly 60 Mgll absorbers having W r (Mgll) > 0.3 A, 
Steidel et al. (1994) identified associated galaxies. Only 16 
of the 45 Mgll systems presented here, of which 21 have 
W r (Mgll) > 0.3 A, have confirmed galaxy counterparts (not 
all fields studied here have been imaged). 

The galaxy properties were obtained from broad-band 
5 (4900/700), 9?(6930/1500), i(8000/1450) and if-band im- 
ages of the QSO fields and their redshifts were spectroscop- 
ically verified to be coincident (within ~ 100 kms -1 ) with 
those of the Mg II absorbers. PSF subtraction of the quasar was 
performed for all fields, enabling small impact parameter galax- 
ies to be identified. Further details of the imaging and spectro- 
scopic observations are described in Steidel et al. (1994). 

The galaxy properties are presented in Table ||; the columns 
from left to right are the quasar field, the galaxy redshift, the 
absolute B magnitude, Mb, the absolute K magnitude, Mk, 
the de-reddened B — K color, the galaxy-quasar impact pa- 
rameter in /i _1 kpc (Ho — 100 km s _1 Mpc -1 , qa = 0.05), and 
a reference if previously published. In a few cases, morphology 
informatio n is available from published HST images (L e Brun 
et al. 1997; |Steidel et al. 1997^ |Steidel 1998[ ). 

It is always possible that a galaxy is misidentified, that an- 
other galaxy or more than one galaxy is giving rise to the ab- 
sorption. Possible selection effects due to misidentifications 
and incompleteness of galaxy redshifts in the individual quasar 
fields are discussed by Steidel et al. (1994) and examined in 
detail by Charlton & Churchill (1996). 

2.2. Data Analysis: Kinematics 

To measure the kinematic "spread" of the low ionization gas 
directly from the flux values, we use the second velocity mo- 
ment of optical depth across the Mgll A2796 profile, defined 

as, 



1/2 



(1) 



where r a (vi) = \n[I c (vi)/ I(vi)] is the apparent optical depth 
meas ured directly from the flux values (S avage & Sembach 
1991) at velocity ~ Av t = (v i - 1 -v i+1 )/2, v t = c(Xi/X obs - 
1), and A fc s = 2796.352(1 + z a b s ). The sums are performed 
only over velocity intervals in which absorption features have 
been detected at the 5 a level, thus eliminating terms consistent 
with noise (see Paper I for details on feature detection). The 
quantity measured by lo v is similar to the standard deviation 
measured from a normal distribution. Th e uncertainty in oj v is 
obt ained from simple error p ropagation (Churchill 1997; also 
see Sembach & Savage 1992). The velocity zero point of each 
system is set at the median wavelength of the apparent optical 
depth distribution of the Mg II A2796 profile. 

The value of lo v depends upon the velocity squared, and thus 
its value is sensitive to the presence of very weak Mgll ab- 
sorption at large velocity. Therefore, u> v is sensitive to the 
equivalent width detection threshold of the HIRES spectra. We 
are 93% complete to a 5 a rest-frame equivalent width detec- 
tion threshold of 0.03 A for the full sample of 45 systems (see 

6 Monte Carlo realizations of our sample could shed light on this issue if, and 
being studied. 



Figure 3 from Paper I). To enforce uniform evaluation of the 
kinematic spread for the full sample, we omitted features with 
^(Mgll) < 0.03 A from the computation of uj v . Six systems 
required this censorship step; they are footnoted in Tablejl]. The 
maximum change in u> v from its uncensored value was a —20% 
difference, except for one of the absorbers (PKS 0454 + 039 at 
z = 1.1532), which was changed by a —40% difference. 

To measure the kinematic "composition" of the low ioniza- 
tion gas, we use Voigt pro file fitting. We use our own program, 
MINFIT (Churchill 1997), which performs a \ 2 minimiz ation 
using the NETLIB-SLATEC routine DNLSl ( More 1978] ). In 
fully saturated profiles, the nu mber of clouds is often under- 
estimated by a factor of a few ( |Churchill 1997 ). A corollary is 
that the cloud column densities, b parameters, and velocities are 
correlative with the number of clouds. We return to this issue 



in S X.2 



3. CHARACTERIZING Mgll ABSORBER PROPERTIES 

We have constructed a subsample of 30 systems, Sample 
CA (CA = Cluster Analysis), in order to examine variations 
and possible trends between the absorption properties listed in 
Table |l| All Sample CA members have measured W r (hya), 
W^Mgll), and uj v and have measurements of or limits on Fell 
and Civ. Systems with upper limits on Vt^Fell) were in- 
cluded in the sample because these limits are already capable 
of demonstrating a paucity of Fell absorption (only seven Sam- 
ple CA systems have Fell limits). Small W r (M.gll) systems 
with upper limits on Civ were excluded from the sample be- 
cause the limits were not always stringent. However, multiple 
cloud systems with limits on W r (Clv) and larger W r (Mgll) 
are included. The members of Sample CA are given in column 
nine of Table fil. 

In Figure 111 we present three dimensional plots of (a) 
W r (Mgll) vs. W r (Lya) vs. W r (Clv), (b) V^.(Mgll) vs. 
W r Qjya) vs. W r (Feii), and (c) W r (Mgii) vs. lu v vs. W r (Civ) 
for Sample CA. Note that distribution of equivalent widths are 
not random; there is a clear trend for W r (Mgll) to increase 
with W r (hya), for example. On the other hand, it is clear 
that there are significant spreads, or variations in the absorp- 
tion strengths. In the case of W^Fell), we see a very large 
range of values (Figure [j]b) that trace W r (Lya), and to a lesser 
extent W r (Mgll). W r (Clv) exhibits a significant spread for a 
given W r (Mgll)-W r (Lyo;) locus (Figure [ija). Moreover, note 
the groupings of W r (Mgll), oj v , and W r (Crv) (Figure |]c). 

From an empirical point of view, the qualitative appearance 
of groupings in Figure [j] would suggest that Mg Il-selected sys- 
tems can be further categorized (quantitatively) by their absorp- 
tion properties. However, our sample is small, having only 
30 data points, and it is not clear whether the apparent group- 
ings would statistically be present in other "realizations" of the 
data 6 . This concern is best addressed by using a multivariate 
analysis in which all available absorption properties are incor- 
porated simultaneously (multidimensional version of Figure [j] 
with one dimension for each property). 

When the Ly a, Fe II, and C IV and Mg II kinematics are si- 
multanously considered, do we find various "classes" of Mg II 
absorbers? In other words, can we quantitatively describe both 
the variations and the trends by considering all the absorption 
line data presented in Table [j]? To address these questions, we 
applied "Tree Clustering" and "if-means Clustering" analysis 

only if, we had a priori knowledge of the distribution functions for all the properties 
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3.1. A Multivariate Analysis 



Multivariate clustering analysis algorithms are designed to 
organize data of many variables into catagories so that natural 
groupings of the data can be examined in a completely unbi- 
ased manner (i.e. no model is imposed upon the data). Full de- 
tails of these techniques can be found elsewhere (e.g. Johnson 
& Wichern 1992; Babu & Feigelson 1996); below, we provide 
limited background material. We used the STATISTICA soft- 
ware package (www.statsoft.com). 

3.1.1. Tree Cluster Analysis 

In tree cluster analysis, the data occupy a multi-dimensional 
space, one dimension for each measured absorption property. 
Clustering algorithms compute the "distances" between each 
pair of points (absorbers) and then amalgamate them into clus- 
ters. We note that tree cluster analysis is not subject to signif- 
icance testing because the result itself is the most significant 
solution under the assumption of no a priori hypothesis regard- 
ing the data. 

The amalgamation process begins with each absorber in a 
unique class by itself and then proceeds by relaxing the cri- 
terion of uniqueness in subsequent steps. With each step, the 
algorithm amalgamates larger and larger clusters of increas- 
ingly disimilar properties, until in the final step all absorbers 
are joined together in a single class. Graphically, clusters ap- 
pear as distinct "branches" in a hierarchical tree, with sim- 
ilarities linked at nodes. We used Euclidean distances [i.e. 
(Eifci - 2/*] 2 ) 1/2 ] and Ward's method (Ward 1963) for the 
amalgamation algorithm. Ward's method minimizes the sum 
of squares of the distances between clusters. The combination 
of Euclidean distances and Ward's amalgamation rule uses an 
analysis of variances approach for which iV(0, 1) standardiza- 
tion of the data is optimal. With this standardization, each vari- 
able is baselined and scaled to have a zero mean and unity stan- 
dard deviation. iV(0, 1) standardization ensures that the dis- 
tance between any pair of points in multidimensional space is 
not biased by a large dynamic range in one or more of the vari- 
ables (dimensions). 



3.1.2. K-means Clustering 

In if-means clustering analysis, we begin with the assump- 
tion of a set number of clusters and the algorithm finds the most 
significant, or distinct, clusters possible. Starting with K ran- 
dom clusters, the algorithm moves points (absorbers) between 
clusters until both the variability within clusters is minimized 
and the variability between clusters is maximized. We find that 
K = 5 is the highest number of clusters allowed such that all 
clusters are considered sign ificant, bas ed upon MANOVA test s 
( lohnson & Wichern 1992 ; also see Mukherjee et al. 1998 ). 
For five clusters, the MANOVA probabilities, p, for accepting 
a result as valid, (i.e. representative of a unique population) are 
highly significant. We obtained p < 1CP 5 for each cluster. 



3.2. Results 

In Figure |[ we present a dendrogram, the tree diagram show- 
ing the results of our cluster analysis. Along the bottom hori- 
zontal axis are the individual Mg II systems identified by quasar 
and absorber redshift. The vertical axis is the linkage distance, 
LD, the distance between clusters. The larger the value of LD 



between two branches, the less related are the objects on each 
branch. 

At LD ~ 16 there is a natural grouping into three clusters. 
These three clusters are predominantly distinguished by very 
strong Lya [right branch, labeled "DLA/Hl-rich"], weaker 
ClV [center branch, labeled "Civ-weak"], and intermediate to 
stronger Civ [left branch, labeled "Civ-strong"] absorption. 

Further insight into the physical differences in the absorbers 
can be obtained if we adopt a "less significant" linkage dis- 
tance. At LD ~ 4 there is a secondary grouping into six clus- 
ters. The Civ-strong and Civ-weak clusters have each been 
separated by variations in Mg II strengths and kinematics and 
the DLA/Hl-Rich cluster has been separated by a spread in the 
Ly a strengths. However, we treat the DLA/Hl-Rich systems as 
a single cluster because the MANOVA significance test from a 
if -means clustering analysis indicates that the further splitting 
of the DLA/Hl-Rich cluster is not justified. On Figure ||, we 
have labeled each of these five clusters as "Classic", "Double", 
"Single/Weak", "Civ deficient", and "DLA/Hl-Rich". 

It is in a diagram of the if -cluster means, Figure ||, that the 
properties distinguishing each cluster become apparent. Across 
the horizontal axis of Figure || are the absorption properties; 
from left to right they are the Mg II, Ly a, C IV, and Fe II absorp- 
tion "strengths", and the "strength" of the kinematic spread. 
Recall that the data have been N(Q, 1) standardized. Thus, for 
each given absorption property, the mean value for all sys- 
tems is zero and the vertical axis is in units of standard de- 
viations. Each cluster is represented by a different data point 
type: Classic (solid circle); Double (solid square); Single/Weak 
(open circle); Civ deficient (solid triangle); and DLA/Hl-Rich 
(solid pentagon). These classes are also connected by unique 
line types (i.e. dash-dot, solid, etc.); it is important to not only 
compare a given property across clusters, but to compare how 
all properties are segregated by cluster. We now briefly discuss 
some distinguishing features of each absorber class. 

3.2. 1 . Classic Systems 

This class appears to have what might be thought of as "typi- 
cal", or non-extreme properties. As seen in Figure |[ this class 
is characterized by having Mg II, Ly a, Fe II, and C IV equivalent 
widths, and a kinematic spread within 0.5 a of the respective 
means for the overall sample. 

We call these systems "Classic" because they can be thought 
to represent the most common type of Mg II absorber observ- 
able in earlie r generation, low resolution surveys (e.g. S teidel & 
Sargent 1992; |Lanzetta, Turnshek, & Wolfe 1987| ). To the sen- 
sitivities of these surveys, it was found that virtually all Mg II 
absorption-selected systems also had Civ absorption. 

3.2.2. ClV "Deficient" Systems 

The C iv-deficient systems have Mg II, Ly a, and Fe II prop- 
erties identical to the Classics and constitute a similar fraction 
of the overall Mg II absorber population. Their distinguishing 
properties are significantly lower C IV absorption strengths. 

3.2.3. Double Systems 

The Doubles are set apart from the Classics by having at 
least twice the Mgll and ClV absorption strengths and Mgll 
kinematics. Since the mean Mgll strength for Doubles is in 
the same regime as of the DLA/Hl-Rich systems, it is the Civ 
strengths and extreme Mg II kinematics that set this class apart 
from the others. Doubles also have Lya and Fell absorption 
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spite of the fact that this Mg II equivalent width is well below 
1.3 A. In other words, even if the Mgll equivalent widths of 
the DLA/Hl-Rich systems were quite small, the class would 
be unchanged in the cluster analysis due to the very strong Hi 
and Fell strengths. This is clearly shown in Figure |[ 

To investigate whether membership in a given class (espe- 
cially the Single/Weak class) may be sensitive to the inclu- 
sion of systems for which only an upper limit is available for 
W r (Clv), we ran the cluster analysis including the four ad- 
ditional systems with information on Lya (without regard to 
whether Civ was measured or an upper limit) A. In this analy- 
sis, the Civ upper limits were treated as measured values. The 
class memberships were unchanged for the original 30 systems. 
The added systems were classified as Single/Weak based upon 
their small W r (Mgll) and narrow kinematics. This was ex- 
pected, given the already large spread in W r (Crv) for the Sin- 
gle/Weak class. 

We thus conclude that the clustering results are robust; they 
are not affected by biases nor are they sensitive to changes in 
the sample. 

3.4. Absorption Strengths 

In Figure |], we present the rest-frame equivalent widths 
(taken from Table [y and from Tables 3 and 4 of Paper I) vs. 
W r (Mgll). We have also included W r (Call) vs. W r (Mgll). 
The panels are ordered by increasing ionization potential from 
the upper left to the lower right and the data point types are the 
same as in Figure |[ 

For the non-DLA/H I-Rich systems, note that the absorption 
strengths of the low ionization species Mgl, Fell, Sill, Alii, 
and Cll increase with increasing W r (Mgll), indicating these 
species arise in the phase giving rise to the Mgll absorption. 
Overall, the higher ionization species Si IV and Civ exhibit 
more scatter, as quantified by the dispersion per observed range. 

We caution that the classification scheme introduced above 
should not be taken to suggest that Mg II absorbers group into 
discretized classes. Discretization is a byproduct of clustering 
analysis. In fact, the distribution functions of the equivalent 
widths plotted in Figure |] are characterized by single modes 
and decreasing tails 7 As such, any single absorption property, 
viewed in this univariate fashion, is distributed continuously. 
However, from the perspective of a multivariate analysis, it is 
clear that the overall properties of Mgll absorbers group in 
well defined regions of a "multi-dimensional space" . 

4. INFERRING IONIZATION CONDITIONS AND KINEMATICS 

To examine the kinematic and ionization conditions, we de- 
fine a new subsample, IC (IC = Ionization Conditions). Sam- 
ple IC includes only those systems with (1) an Mgll equivalent 
width detection threshold less than or equal to 0.02 A, and (2) 
no unresolved saturation in at least four adjacent pixels (1.33 
resolution elements) in the Mgll A2796 profile. These selec- 
tion criteria enforce a high level of accuracy in the number of 
Voigt profile components, N c i, and their column densities, ve- 
locities, and b parameters. Simulations of blended, multiple 
component Mgll profiles with these characteristics show that 
the distribution of Voigt profile parameters output from Voigt 
profile fitting is consi stent (99% confi dence) with those used to 
generate the profiles (Churchill 1997). Data with lower signal- 
to-noise ratios result in a slight paucity of Voigt profile com- 

7 The exceptions are Wr(Mgl), W r (Fell), and W r (Lyot), which are bimodal due to the DLA/Hl-Rich class. However, it is not clear if this bimodality is due to 
small numbers and/or is due to the relatively large number of DLAs in our sample. 



strengths that are systematically greater than those of the Clas- 
sics. 

3.2.4. Single/Weak Systems 

The Single/Weak class is defined foremost by very small 
Mgll strengths and kinematics. In fact, with the exception of 
two systems, all are single, narrow clouds (often unresolved in 
the HIRES spectra). 

The Single/Weak systems are underrepresented in Sample 
CA, mostly because seven of the systems in our database have 
no measured Lya. If these seven systems were included with 
their Civ limits artificially treated as detections, the mean Civ 
for this class would drop to —1.1 (this being an upper limit) 
in the K-means cluster diagram (Figure ||). Thus, there is a 
spread in the distribution of C IV strengths associated with Sin- 
gle/Weak systems. The true mean Fell is also probably lower 
than that shown on Figure || because the limits on Fe II in the 
HIRES spectra scatter about the measured mean. Thus, higher 
sensitivity spectra would reduce the mean Fell for this class 
and could reveal a spread in Fell strengths. Only for the Sin- 
gle/Weak class, for which many of the properties were below 
our detection thresholds, are the above issues pertinent. 

3.2.5. DLAMl-Rich Systems 

The DLA/Hl-Rich class has very strong Mgll, and ex- 
tremely strong Lya and Fell, relative to the overall Mgll ab- 
sorber population. They have C IV strengths below that of the 
Classics and Mg II kinematics typical of the Classics and C IV- 
deficient absorbers. Their Mg I and Fe II strengths are five to ten 
times greater than those of the typical Classic absorber. Unlike 
the elevated Lya strengths, which are due to broad damping 
wings, the Mg I and Fe II strengths are driven by the line-of- 
sight kinematic spreads of the gas for multiple, saturated com- 
ponents. For the intermediate ionization species, Sill, Alii, Cll 
and Si III, the strengths are somewhat greater than those of the 
Classics, on average. The higher ionization species, Si IV, N V 
and O VI (to the extent the latter two have been measured), have 
strengths consistent with those of the Classics absorbers. 

3.3. Robustness of the Clustering Analysis 

Two concerns regarding the clustering analysis are: (1) pos- 
sible biasing due to an slight, but apparent, overabundance of 
absorbers with W r (Mgll) > 1.3 A, and (2) the exclusion of the 
"weak" systems for which only upper limits on the C IV equiv- 
alent widths were measured. 

The systems with W r (Mgll) > 1.3 A are all members of 
the DLA/Hl-Rich class. Of these, four are bonified DLAs with 
AT (Hi) > 2 x 10 20 cm" 2 , whereas two have column densities 
below this classical threshold (i.e. they are Hi-Rich systems). 
In an unbiased survey for z < 1.65 DLAs, using Mgll absorp- 
tion as a selection method, Rao & Turnshek (1999) found that 
14% of the Mgll systems with W r (Mgll) > 0.3 A are DLAs. 
In our sample, 22 systems have W r (Mgll) > 0.3 A; thus, 18% 
of our sample is comprised of DLAs. This is not inconsistent 
with the unbiased results of Rao & Turnshek. 

The presence of a DLA/Hl-Rich class in our sample is not 
due to the slight bias (overabundance) of the larger W(Mgll); 
it is mostly defined by large Hi and Fell equivalent widths. 
In fact, one DLA has W r (Mgll) = 0.9 A (z = 0.5764 to- 
ward Q 0117 + 213), and was classified as DLA/Hl-Rich in 
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ponents, and therefore (to compensate) the resulting compo- 
nents have column densities and b parameters that are too large. 
Profiles with severely saturated cores (e.g. DLA/Hl-Rich ab- 
sorbers) systematically have fewer Voigt profile components by 
a factor of three and have unconstrained column densities and 
b parameters. The equivalent width threshold criterion directly 
translates to a signal-to-noise ratio criterion of S/N = 22 per 
resolution element in the continuum. Sample IC membership is 
given in column nine of Table [j]. 

4. 1 . Ionization Conditions 

4.1.1. Single/Weak Systems 

Normally, it is difficult to interpret the ionization conditions 
based upon equivalent widths because simple curve of growth 
arguments are muddled by the possibility of unresolved satura- 
tion in multiple component absorption. However, these argu- 
ments can hold for the Single/Weak clouds. 

As shown by Churchill et al. (1999a), under the assump- 
tion of photoionization, Single/Weak clouds with W^Mgll) < 
0.15 A, cannot give rise to W r {C\w) > 0.2 A (see their Fig- 
ure 12). The upper limit on W,.(Civ) is significantly more 
restrictive when Fell is detected in the Mgll cloud (i.e. the 
cloud is constrained to have low ionization conditions). Based 
upon this analysis, we find that roughly half of the observed 
Single/Weak absorbers likely have multiple ionization phases, 
with Civ and some portion of the Lya absorption arising in 
spatially distinct, higher ionization material. These findings 
are c onfirmed in a thorough study using CLOUDY (F erland 
1996) photoionization models tuned to the measured Mgll col- 
umn densities and constrained by the full complement of ion- 
ization sp ecies and transitions covered in FOS spectra (R igby 
et al. 1999) 

4.1.2. Multicloud Systems 

We now use photoionization modeling to investigate whether 
multiphase ionization can also be inferred for the multicloud 
Classic, C iv-deficient, and Double systems. We consider 
whether the measured C IV and/or Si IV absorption can or can- 
not all arise in the Mg II clouds, even under extreme assump- 
tions about their ionization conditions. The following is based 
upon the methods employed by Churchill & Charlton (1999) in 
their study of the z = 0.9276 systems toward PG 1206 + 459. 

We assumed that the Mgll clouds in a given system are 
in p hotoionization equilibrium and used CLOUDY (F erland 
1996) to model their ionization conditions. A Haardt & Madau 
(1996) extragalactic spectrum at z — 1 was used for the ioniz- 
ing radiation. Based upon experiments with various single and 
extended starbu rst galaxy spectral energy distributions (B ruzual 
& Chariot 1993) with assumed high photon escape fractions, 
we determined that CLOUDY models involving Mgll, Fell, 
Si IV, and C IV are not highly sensitive to the chosen spectrum 
for z ~ 1 ( |Churchill & Charlton 1999| ); thus our general con- 
clusions are not sensitive to the assumed spectral energy distri- 
bution. 

We obtained a predicted C IV and Si IV equivalent width for 
each system by tuning a CLOUDY model to the observed Mg II 
and Fell column densities in each of its clouds. We then mea- 
sured the absorption strengths in synthesized HIRES spectra us- 
ing the velocities and b parameters thermally scaled from those 
measured for the Mg II clouds [see Churchill & Charlton (1999) 
for details.] 

Since our goal is to infer if multiphase ionization structure 



is possibly present in a given system, we have forced the mod- 
els to yield the maximum amount of Civ and Si IV that can 
arise in the Mgll clouds. If a cloud has detected Fell, then, 
for a given [a/Fe] abundance pattern under the assumption of 
negligible dust depletion, the Fell and Mgll column densities 
uniquely determine the ionization parameter (the logarithm of 
the ratio of the number density of hydrogen ionizing photons to 
the number density of hydrogen), nearly independent of metal- 
licity. The ionization parameter uniquely dictates the C IV and 
Si IV column densities in the clouds. Typical ionization param- 
eters in the clouds with Fell are —4 to —3. We use a solar 
[a/Fe] abundance pattern, since this yields a higher ionization 
parameter for a given ratio of Mg II to Fe II column densities, 
and thus yields a maximized Civ and Si IV strength. If there is 
no detected Fe II for a cloud (either a limit or no coverage), then 
the ionization parameter is pushed to as high a value as possible 
without the cloud size exceeding 10 kpc. Typically, this occurs 
at an ionization parameter of ~ —1.9. A solar metallicity was 
assumed for all clouds, again in order to maximize the predicted 
Civ and SilV strengths (for a high metallicity cloud, the hydro- 
gen column density required to give rise to the observed Mg II is 
smaller, and the ionization parameter can be pushed to a higher 
value). We did not apply any other constraints from the data 
(e.g. Sill, Cll, etc.) than those described above. 

In Figure ^, we present the predicted maximum equivalent 
widths vs. the observed equivalent widths for Civ (left panel) 
and Si IV (right panel). The mean error in the measured val- 
ues is depicted by the open-box data point in the upper left 
of each panel. Diagonal lines demarcate those systems that 
are inferred to have multiphase ionization (lower right) from 
those that could be single phase. The maximum predicted Civ 
strength that can arise in photoionized Mg II clouds is more than 
3 a short of the observed W^Clv) for a number of systems. 
Even two of the C iv-deficient systems are below the line. This 
implies that even though these particular systems have below 
average C IV absorption strengths, the C IV may still arise in a 
distinct, higher ionization phase and not in the Mgll clouds. 
For Si IV, it appears that fewer systems would require multi- 
phase structure in this ionization species; Si IV likely resides in 
the Mg II clouds in a larger fraction of Mg II selected absorbers 
than does the Civ. 

4.2. Mg II and C IV Kinematics 

4.2. 1 . Number of Clouds 

In Figure |J we show these same absorption strengths vs. the 
number of clouds, or subcomponents, N c i, obtained from Voigt 
profile decomposition. The data point types are the same as for 
Figure fj| 

There is a strong correlation of W r (Mg II) with N c i, the num- 
ber o f Voigt profile components (also see P etitjean & Bergeron 
1990). A linear least-squares fit, which we have drawn as a dot- 
ted line in Figure U yielded a slope of 0.069 ±0.005 A cloud" 1 . 
Note that this relationship is likely resolution dependent. Also 
note that the fit is strongly influenced by the two Double sys- 
tems with large N c i. 

We find that there is not as "linear" a dependence for 
W^Fell) nor for W r (Clv) with Mgll kinematics. There is, 
however, a tight correlation of W r (Lya) with the number of 
Mgll clouds. We obtained a linear least-squares fit (drawn as 
a dotted line) with slope 0.57 ± 0.08 A cloud -1 . The linear 
relationship for W^Mgll) vs. N c i suggests that majority of 
neutral hydrogen equivalent width arises in the Mgll clouds 
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themselves. 



4.2.2. Mg II Kinematics 



In Figure |], the absorption strengths of Mgll, Fell, Lya, and 
Civ are plotted vs. the Mgll kinematic spread, lo v . The data 
point types are the same as for Figure [£| 

In Figure |J, the scatter in W r (Mg II) arises because the Mg II 
equivalent width is dominated by the largest clouds, which are 
usually clustered in subgroups with velocity spreads less than 
- 20 km s" 1 [see Figure 12 of Charlton & Churchill (1998)]. 
The kinematic spread, on the other hand, is sensitive to the pres- 
ence of small W r (Mg II) clouds with large velocities (u> v oc v 2 ). 
The "limit" along the upper left of the data, shown as a dotted 
line, is due to saturation in the profile line cores; there is a max- 
imum equivalent width for a given velocity spread when the 
profile is black bottomed. The "limit" along the lower right is 
due to detection sensitivity; for the signal-to-noise ratios of the 
sample, there is a minimum detectable W r (Mgll) for a given 
velocity spread. Higher quality data would be required to de- 
termine if there is an actual paucity of systems with large kine- 
matic spreads and extremely weak Mg II absorption. 

The strong correlation between W r (Crv) and u v has previ- 
ously been discussed by Churchill et al. (1999b). The corre- 
lation is driven by the three Double systems (z = 0.8514 to- 
ward Q 0002 + 051, z = 0.9110 toward PKS 0823 - 223, and 
z = 0.9276 toward PG 1206 + 459) and the two kinematically 
extreme Classic systems (z — 1.3250 toward PG 0117 + 213 
and z = 0.7908 toward PKS 2145 + 067), which all have uo v > 
75 kms -1 . We have drawn in the maximum-likelihood lin- 
ear fit to the data, which has slope 15.4 ± 0.2 mA (km s -1 ) -1 . 
More data on Double systems at intermediate redshifts, selected 
by W r (MgIl) > 0.6 A, would be useful for determining how 
tight the relationship between W r (Clv) and Mgll kinematics 
remains for the largest kinematic spreads. 

Three of the five Classic systems above the correlation line 
are from the lowest signal-to-noise spectra; there could be 
missed, small equivalent width components resulting in a small 
uo v . Even so, there is a significant scatter in W r (C\v) at a 
given uo v for u> v < 60 km s _1 . The Single/Weak systems have 
M^.(Clv) ~ 0.5 A and uj v ~ 5 kms -1 , whereas the ClV- 
deficient systems have W r (Civ) < 0.2 A and u> v ~45kms _1 . 

4.2.3. ClV Kinematics 

In Figure || we display the suite of C IV doublets feasured in 
FOS spectra (we exclude ground based data; see Paper I) that 
are not blended with other transitions. The velocity window for 
the Civ doublets is 2000 km s _1 . Above each doublet member 
are ticks giving the velocities of the Mg II Voigt profile compo- 
nents. The Mgll A2796 profiles are also shown with the Voigt 
profile models superimposed on the data, which are ordered by 
increasing kinematic spread, u v . The velocity window for the 
Mgll data is 600 kms -1 . 

To demonstrate which profiles are resolved, and to quan- 
tify the degree to which they are resolved, we fit the indi- 
vidual Civ members with single Gaussians while holding the 
width constant at the value of the instrumental spread function, 
a = 98 km s~ x . We have superimposed these unresolved fits 
on the Civ AA1548, 1550 profiles. We note that several of the 
spectra were obtained prior to the refurbishing mission in which 
COSTAR was installed. For "large" aperture aquisition modes 

8 Interestingly, this DLA has W r (Crv) = 1.63 A, twice that of the DLA/Hl- 
Sargent 1992) 



(1" slit), the instrumental spread function of the pre-COSTAR 
instrument has extended wings. PKS 0454 — 220 is the only 
object acquired in this mode; thus, the unresolved fit to the 
z — 0.4744 system should be viewed with some discretion. 
Several of the unresolved fits are consistent with the C IV data, 
suggesting that the adopted instrumental spread function accu- 
rately represents the pre-COSTAR and post-COSTAR instru- 
ment for the smaller apertures modes. 
We computed the quantity, 



F r = 



W (measured) — ^(unresolved) 
^(measured) 



(2) 



where ^(measured) is the measured Civ equivalent width 
adopted for Table [I], and ^(unresolved) is the Civ equivalent 
width from the unresolved fits. Thus, F r measures the frac- 
tion of the integrated flux due to resolved velocity structure in 
the C IV doublet. We have noted the F r value and its 1 a un- 
certainty for each Civ profile in Figure [| There is a 3 a cor- 
relation between F r and uj v . We have thus demonstrated that 
the tight W r (Crv) vs. oj v correlation is due to resolved velocity 
structure in the C IV absorbing gas. 

Velocity structure is not unexpected. Observed at higher res- 
olution, C IV profiles from Milky Way halo gas exhibit veloc- 
ity structure over a ~ 200 km s _1 spr ead, though the compo - 
nents are typically blended together ( |Sembach et al. 1999b ). 
For higher redshift galaxies, the C IV profiles als o show distinct 
velocity structure (Petitjean & Bergeron 1994) with velocity 



spreads up to 200 km s _i . In Figure |?| we present an exam- 
ple of a kinematically complex Civ doublet (~ 9 km s _1 reso- 
lution) and two Fe II transitions (Mg II was not covered) for the 
damped Ly a absorber at z — 1.7763 toward Q1331+170 8 . The 
ticks above the C IV profiles were obtained by a crude Voigt pro- 
file fit to the strongest components to parameterize the kinemat- 
ics; we used this fit to synthesize the Civ doublet as it would 
be observed with FOS resolution and pixel sampling (infinite 
signal-to-noise ratio). This is shown in the separate, bottom 
panel of Figure^, with the Civ A1548 and A1551 component 
ticks above the continuum (not the Mgll ticks as in Figure |8|). 
An unresolved fit is superimposed on the synthetic FOS doublet 
for which we measured F r = 0.27 ± 0.02. This value is com- 
parable to F r for the z = 0.5505 system toward Q 1241 + 174 
and the z = 0.9110 double system toward PG 0823 - 223. 

In cases where the Mgll profiles are highly asymmetric, it 
appears that the C IV profiles exhibit similar asymmetry in the 
same sense. That is, the flux discrepancies between the unre- 
solved fits and the C IV profiles are aligned in velocity with the 
kinematic outliers seen in Mg II. Since we have demonstrated 
that, in most cases, the C IV does not arise in the same phase as 
the Mgll clouds (see § 4.1.2), we infer that the ClV often arises 



in a physically distinct structure from the Mg II, but is aligned 
kinematically with the Mg II. 

5. ABSORPTION-GALAXY PROPERTIES 

We tested for correlations between the galaxy properties and 
all the absorption line data presented in Paper I, including 
all possible combinations of equivalent width ratios. We per- 
formed Kendall and Spearman non-parametric rank correlation 
tests using the program ASURV (L a Valley, Isobe, & Feigelson 
1992). The ASURV algorithm provides for measurements that 
are either upper or lower limits. 

Rich systems in our lower redshift sample. It has Wr(Mgll) = 1.3 A (Steidel & 
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Out of 114 tests performed, seven resulted in correlations 
above the 2.5 a level 9 ; these are: W r (MgIl), Mgll dou- 
blet ratio, and W r (Silv)/W / r (Clv) vs. impact parameter; 
W r (Cll) and W r (Mgll) vs. M K ; and W r (Ly a)/W r (Clv) and 
W r (Sill)/W r (Sim) vs. B — K. When DLA/Hl-Rich systems 
were removed from the sample, the significance levels dropped 
slightly below 2.0 a. No test resulted in a correlation at a 3.0 a 
or greater significance level. 

5.1. Global Ionization and Density Structure ? 

In Figure [l(] we present plots of the absorption properties 
used in our multivariate analysis, u> v , W^Mgll), W r (Feu), 
W r (Lya), and Vt^Crv), vs. host galaxy properties. 

W r (Mgll) is anti-correlated with impact parameter (de- 
creases with galactocentric distance) with a 2.7 a significance 
and is correlated with Mk (increases with decreasing luminos- 
ity), also with a 2.7 a significance. Both these trends arise be- 
cause DLAs, which give rise to the largest W r (Mgll), are often 
observed at low impact parameters and often have low lumi- 
nosities (M K > -24). 

The W^Crv) vs. K luminosity relationship can provide a 
test for halo models in which a virialized hot phase (T ~ 10 6 K) 
pressure con fines the low ionization phase (e.g. M o & Miralda- 
Escude 1996). These models predict an anticorrelation between 
W r (Cw) and galaxy mass (i.e. K luminosity); C IV is predicted 
to be enhanced in smaller, lower mass galaxies, which have low 
gas pressures in their halos. As seen in Figure |IoJ there is a vi- 
sual, yet not statistically significant, trend for the most massive 
galaxies to have smaller W r (Clv) in our small sample. This is 
driven by the high luminosity [and very red; see W r (Crv) vs. 
B — K] C iv-deficient systems. 

In Figure [H], we have plotted the equivalent width ratios of 
Si iv to Civ, Lya to Civ, Lya to Mgll, Fen to Civ, and 
the Civ and Mgll doublet ratios vs. impact parameter. We 
have selected combinations only of the gas properties used 
in our multivariate analysis, with the exception of the ratio 
W / r (Silv)/W / r (Clv). This ratio is anticorrelated with impact 
parameter at the 2.8 a level when all galaxies are included and 
at the 2.1 er level when DLA/Hl-Rich systems are excluded. If 
this trend were to hold for a larger sample, it would be consis- 
tent with the scenario in which halo gas is photoionized by the 
extragalactic UV background flux resulting in a increasing ion- 
ization level with increas ing galactocentric distance (see Sav- 
age, Sembach, & Lu 1997). 

The Mg II doublet ratio is correlated with impact parameter 
at the 2.7 a significance level. This reflects the fact that the 
gas tends to become optically thin at larger galactocentric dis- 
tances, where W r (Mgll) is systematically smaller. Again, the 
DLAs dominate the trend. 

The ratio W r (Lya) /W r (Mg II) is not correlated with impact 
parameter when all galaxies in the sample are included. How- 
ever, when the DLA/Hl-Rich systems are excluded a 2.5 a 
level correlation is observed. This is governed by the gen- 
eral decrease in Mg II equivalent widths with impact parame- 
ter, since Lya shows no trend with impact parameter (when 
DLA/Hl-Rich systems are excluded). This particular trend 
is difficult to interpret because of the very different curve of 
growth behavior for Mg II and Ly a. It is not possible to address 
the possibility of a metallicity or ionization gradient. The loca- 
tion of the Lyman limit break was covered for only six of the 
galaxies and was detected for all six, sampling impact parame- 



ters from 10 to 40 kpc. 

5.2. Ionization Conditions and Colors 

In Figure |l2[ we have plotted the C IV to Mg II, Ly a to C IV, 
and Lya to Mg II equivalent width ratios vs. host galaxy B — K 
color. Both W r (Lya)/W r (Clv) and W r (Lya)/W / r (Mgll) are 
correlated with B — K at a 2.6 a significance level. This trend 
in our small sample is governed by both the DLA/Hl-Rich sys- 
tems and the C iv-deficient systems being associated with red- 
der galaxies. 

5.3. Correlation Test: Overall Results 

The degree of scatter in the absorption properties presented 
here are similar to that found by Churchill et al. (1996) for the 
Mg II absorption properties only. None of the presented absorp- 
tion properties correlated with the galaxy properties at a high 
significance level. This is not to say that definite trends, or even 
statistically significant correlations, do not exist. The large scat- 
ter in the properties signifies large local variations from line of 
sight to line of sight through the galaxies; however, even with 
our small sample, we have uncovered some suggestive global 
trends. 

We also find that the inclusion or exclusion of DLA/H I-Rich 
systems in any given test significantly alters the statistics. The 
host galaxies of DLAs are seen to have a wide range of lumi- 
nosities, morphologies, colors, and surface brightnesses (e.g. 
Le Brun et al. 1997 ; Rao & Turnshek 1998 ), and do not repre- 
sent a homogenous population, despite the fairly homogeneous 
absorption properties by which they are selected. 

6. DISCUSSION 

We reemphasize that a taxonomic scheme based upon a mul- 
tivariate clustering analysis will, by its very nature, result in a 
discretization of what may actually be a continuum of proper- 
ties. Nonetheless, if w e adopt the clustering analysis results at 
face value (see § 3.2) and segregate the wide variety of Mgll 



absorber properties into individual "classes", we may be able 
to find order in an otherwise complex array of gas properties 
and kinematics. The hope is to gain further insights into defin- 
ing distinctive or "characteristic" properties of Mg II absorbers 
and to understanding the galactic processes that give rise to the 
observed range of these "characteristics". 

What follows is primarily a speculative discussion, given the 
fact that there is still very little data upon which inferrences can 
be based. 

6.1. Classic vs. C TV-deficient Systems 

The existence of C iv-deficient systems with W r (Clv) < 
0.15 A, which have kinematics similar to the Classics, indicates 
that the connection between C IV absorption and Mg II kinemat- 
ics does not operate the same in all galaxies (or at all locations 
in galaxies). 

There are only three C iv-deficient systems with measured 
galaxy properties {z = 0.4297 toward PKS 2128 - 123, 
z = 0.6601 toward Q 1317 + 277, and z = 0.7291 toward 
PG 0117 + 213). The host galaxies of these three systems 
are reddish (B — K = 3.3, 3.8, and 4.0) and are among the 
more luminous (Mk = —24.7, —25.7, and —26.3). Also, they 
are probed at large impact parameters (D — 32, 38, and 36 
hr 1 kpc), respectively (see Figure [l(J). The redder colors and 
large luminosities would suggest massive, early-type galaxies. 



3 For 1 14 correlations tests, only one test should result in a 2.5 a or greater significance level at random. 
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Based upo n HST image s of the z = 0.6601 and z = 0.7291 
absorbers (Steidel 1998), we roughly (visually) classify them 
to be an edge-on SO galaxy (probed on its minor axis) and a 
face-on SBa galaxy, respectively. 

Despite small numbers, these relatively high masses, early- 
type morphologies, and large impact parameters may be a clue 
to the observed spread in W r (C IV ) for a given Mg II absorp- 
tion profile type. Semianalytical models of pressure confined 
gaseous halo s predict smaller Civ strengt hs in more massive 
galaxies (e.g. Vlo & Miralda-Escude 1996), and such a trend is 
not inconsistent with our data. The weak C IV cannot simply be 
due to the line of sight through the galaxy, because some Clas- 
sics at large impact parameters have very strong Civ strengths. 
Also, based upon cross-sectional a rguments fro m the number 
of Civ absorbers per unit redshift ( |Steidel 1993 ) and upon the 
expected lowe r pressures and gas densities a t larger galactocen- 
tric distances ( Mo & Miralda-Escude 1996) ), it is expected that 
Civ absorption would be strong at large impact parameters. 

Perhaps strong C IV is expected in regions of relatively pro- 
nounced star formation, which would be consistent with C IV- 
deficient absorbers being associated with earlier-type galax- 
ies. Alternatively, in C iv-deficient absorbers, the Civ could 
be ionized away due to collisional processes, such as with 
the intragroup medium scenario proposed by Mulchaey et al. 
(1996). If so, strong Ovi absorption would be predicted; un- 
fortunately, we cannot address Ovi absorption for the ClV- 
deficient systems [only one is "measured" to have an unre- 
strictive upper limit on W r (Ovi)]. A further possibility is 
that galaxies with extended X-ray emission , usually early-type 
galaxies (e.g. Mathews & Brighenti 1998), may preferentially 
be Civ- and Ovi-deficient because of the extremely hot envi- 
ronment. If so, the Mg II absorbing gas would necessarily arise 
in Hi gas, which is of external origin for ellipticals from galaxy 
mergings ( Knapp, Turner, & Cunniffe 1985; van Gorkom 1992) 
and is sometimes found to be in disk-like structures extending 

5 to 10 t imes the optical radii in "dust-lane" ellip tical galax- 
ies (e.g. Morganti, Oosterloo, & Tsvetanov 1998| ) Thus, one 
might hypothesize that Mg II absorption associated with ellipti- 
cal galaxies arises preferentially in post-merger ellipticals. 

6.2. Classic vs. Double Systems 

There are at least two obvious explanations for the Double 
systems; they might be two Classic systems clustered in line- 
of-sight velocity (bound or unbound), or a primary galaxy and 
a satellite (e.g. possible an interacting LMC-like object). 

The three Double systems in our sample have the largest 
W r (Cl\), which are consistent with that of two typical Clas- 
sic systems. Alternatively, they could arise in a primary galaxy 
either undergoing a minor merger or residing in a group with 
several minor galaxies. Using the Local Group as a model and 
applying the simple cros s-sectional d ependence for Wr(MgIl) 
with galaxy lu minosity ([Steidel 1995t also see McLin, Giroux, 

6 Stocke 1998), we estimated the probability of intercepting a 
"double" absorber for a random line of sight passing through 
the Milky Way (line-of-sight kinematics were not considered). 
We found a 25% chance of intercepting both the LMC and the 
Milky Way, and a 5% chance of intercepting both the SMC and 
the Milky Way. All other galaxies in the Local Group have 
negligible probabilities of being intercepted for a line of sight 
passing within 50 kpc of the Milky Way. 

An alternative scenario is that Double systems are otherwise 
typical Classic systems, but for the interception of kinematic 
outlier clouds, either due to infalling or ejected "high velocity" 



gas. This picture would naturally invoke elevated star forma- 
tion to explain the large W r (Civ) in Doubles and, possibly, the 
correlation between W r (Clv) and Mg II kinematics. Inf alling 
material can enhance star formation (iHummel et al. 199(| Lutz 



1992; |Hernquist & Mihos 1995| ; |Hibbard & van Gorkom 19960 
which can prov ide the energetics to support a multiphase halo 
( |Dahlem 1998 ) with strong Civ absorption [i.e., as with the 
Galaxy (Savage et al. 1997)]. Consistent with this scenario 
would be that the bluest galaxies would preferentially be as- 
sociated with Double systems. Only one Double system has 
measured galaxy properties (z — 0.8519 toward Q 0002 + 051), 
and we note that its color (B — K = 2.9) is among the bluest 
in the s ample, consi stent with a star forming object. An HST 
image (Steidel 1998) reveals the galaxy to be compact with a 
spherical, featureless morphology (no evidence for merging). 

6.3. Classic vs. DLAMl-Rich 

Though the Classic Mg II systems are well u nderstood to se- 
lect normal, bright galaxies ( Steidel et al. 1994 ), the DLA/ sys- 



tems are seen to select [based upon N (Hi) > 2 x 10 cm" ] 
an eclectic population o f low luminosity, lo w surface bright- 
ness, mid^waxf£alaxies^e^^m^^T99^Rao & Turnshek 
1998, 1999). This might suggest that not all classes of Mg II ab- 
sorber have a systematic connection with galaxy type or galaxy 
properties. 

What is interesting, however, is that the Mgll kinematic 
spreads for DLA/H I-Rich systems are tightly clustered around 
~ 60 km s _1 , have no weak, kinematic outlier clouds, and have 
intermediate to weak C IV absorption strengths. This is in con- 
trast to those at redshifts 2-3, where strongly absorbing, higher 
velocity outlying clouds and often "dou ble" profiles ar e charac- 
teristic of low ionization DLA profiles (Lu et al. 1996). 

Whatever physical process or environment gives rise to kine- 
matic outlying Mg II clouds and strong Civ, it is apparently not 
acting in intermediate redshift DLA/H I-Rich systems. Alterna- 
tively, these trends could be due to selection effects of sightlines 
through Hl-rich environments, which implies a line-of-sight 
dependency for absorber class. 

6.4. Absorber Class Evolution? 

We have shown that there is a signicant correlation be- 
tween W r (Cw) and the Mgll kinematic spread. But, exactly 
how do the Mgll kinematics vary with W r (Clv) for a given 
W r (Mgll)? In Figure ||J we have plotted the Mgll A2796 
profiles from our HIRES/Keck spectra in the approximate lo- 
cations they occupy in the W r (Clv)-W r (Mgll) plane. Based 
upon this schematic, we have an impression of how the kine- 
matic "complexity" of the absorption profile relates to both the 
Mgll and Civ absorption strengths. Roughly, each of the five 
classes obtained from our multivariate analysis falls in a well 
defined region, which we have marked to guide the reader. Note 
that, in general, the C IV strength is largest when the kinematic 
"complexity" is the greatest; Civ absorption appears to have 
no direct connection with the optical depth of the strongest 
Mgll component, nor with the total number of Voigt profile 
components. Most remarkably, it would appear that the gross 
characteristics of the Mg II kinematics can be predicted simply 
based upon the absorber's location on the W r (Clv)-W r (Mgll) 
plane. 

Since the distribution of Mg II e quivalent widths evolves with 
redshift (Steidel & Sargent 1992), and this must be reflected in 
the Mg II kinematics as well, it seems reasonable to assume that 
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ing "double" black-bottomed profiles. 

Of the systems shown in Figure |l5j only the z — 1.7945 
system toward B2 1225 + 317 is represented on Figure |l4| . 
Though the Civ strengths of the other systems shown in Fig- 



ure 
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evolution should be discernable on the W r (Clv)— W r (Mgll) 
plane. 

There are at least two possible types of evolution of absorber 
classes: ( 1 ) the number density per unit redshift of a given class 
could evolve, either diminishing with or increasing with time, 
and/or (2) any absorber in a given class could evolve into an- 
other class. The former is systematic and would indicate a 
global, cosmic evolution, which might imply the existence of 
absorber classes not seen in our intermediate redshift sample. 
The latter form of evolution would be stochastically related to 
the processes occuring in galaxies and their environs. 

We present the observational Vt / , (Clv)-H / r (Mgll) plane in 
Figure [l4]. The intermediate redshift data (taken from this 
study) have small, filled circle data points. We have included 
lower and higher redshifts data from the literature, which we 
have listed in Table || The higher redshift absorbers (1.2 < 
z < 2.1) are taken from Bergeron & Boisse (1984), Boisse 
& Bergeron (1985), Lanzetta et al. (1987), and Steidel & Sar- 
gent (1992) [six-pointed stars]. The lower redshift points 
(0.1 < z < 0.6) are taken from Bergeron et al. (1994) [down- 
ward pointing, open triangles], and the z ~ data 10 are 
taken from Bowen et al. (1995, 1996), Jannuzi et al. (1998), 
and Bowen (1999, private communication) [open diamonds]. 
There are four extraordinary systems listed in Table || that have 
W r (Mg II ) ranging from ~ 4 to ~ 7 A. These have not been 
presented in Figure [l4| 

Three points from Figure [l4] are that (1) there are absorbers 
in the higher, lower, and z ~ redshift samples with Civ 
and Mgll strengths typical of each of the five classes found 
in our clustering analysis of the intermediate redshift systems; 
(2) there are both higher and z ~ redshift absorbers occupy- 
ing regions of the W r (C l\)-W r (Mg il) that are not represented 
in our intermediate redshift sample; and (3) as with the inter- 
mediate redshift systems, the higher redshift systems exhibit a 
large range of W r (Clv) for a given W r (Mgll). Though there 
is no suggestion for such a spread in the lower redshift systems, 
there are too few measurements to characterize the spread in 
W r (Clv) values. 

Given the fairly systematic dependence of kinematics in the 
W r (Clv)-W r (Mgll) plane, we can infer that the Mgll kine- 
matics of the lower and higher redshift systems that occupy lo- 
cations consistent with a given absorber class are likely to be 
similar to that class' kinematics. 

The higher redshift systems with W r (MgIl) > 2 A, are a 
class of absorber not present in our sample. Note the large 
spread in W r (C IV ) for these systems, which would suggest that 
some are "C IV deficient". This may be indicative that the phys- 
ical processes giving rise to the range of W r (Clw) in interme- 
diate redshift systems are also occuring in these W r (Mg n) > 
2 A, higher redshift systems. 

In Figure [TJ, we present the HIRES/Keck Mgll A2796 pro- 
files of fourz > 1 Mgll systems with W r (Mgll) > 2.0 A. Note 
that the profiles are optically thick with virtually unity doublet 
ratio over a the full velocity span, which is typically 250 km s _1 
or greater. We also show the z ~ Mg II profile measured in the 
spectrum o f SN 1 993 J in M8 1 using the GHRS on HST ( Bowen 
et al. 1995; spectrum provided courtesy of Dr. D. Bowen). In 
comparison to the smaller equivalent width, intermediate red- 
shift absorbers, the Mgll kinematics and absorption strengths 
of these W r (Mg n) > 2 A systems are clearly unique, exhibit- 

10 We note that the z ~ sample is small and is selected based upon a bright, background quasar having a small projected separation from a nearby galaxy, as 
opposed to being Mgll-absorption selected (Bowen et al. 1995); this may have introduced a bias, possibly related to impact parameter distribution, as compared to 
the higher redshift data. 



are unmeasured or unpublished, it is clear that these sys- 
tems constitute their own unique "class" (with respect to the 
five classes found in our clustering analysis). As shown in Fig- 
ure Ol the Mg II kinematics are suggestive of "double-DLAs", 
and, in fact, the z = 1.7945 system tow ard B2 1225 + 317 is 
a DLA ( [Bechtold, Green, & York 1987| ). In high quality data, 
"Double-DLAs" might be expected to have slightly asymmet- 
ric damping wings, reflecting the different column densities in 
the two systems. 

Steidel & Sargent (1992) found that the number density of 
these large equivalent systems decreases with redshift over the 
range 0.3 < z < 2.2 (the full range over which the Mgll dou- 
blet is observable from the ground). Based upon the SN 1993J 
spectrum of M8 1 and the Galaxy in absorption (the line of sight 
passes through half the disk and halo of M81, half the disk and 
halo of the Galaxy, and through "intergalactic" material appar- 
ently from the strong dwarf-galaxy interactions taking place 
with both galaxies), one possibility is that these systems arise 
from two galaxies with low impact parameters that happen to 
have a line-of-sight superposition. 

For this scenario to be consistent with the Mgll absorber 
evolution, the number of galaxy pairs would need to decrease 
in step with the evolution of the absorbers themselves over 
the same redshift regime. In fact, over the redshift interval 
1 < z < 2, it is seen that the galaxy pair fraction, where 
pairs are defined to have projected separations less than 20 kpc, 
evolves propo rtional to (1 + z) p , with 2 < p < 4 ( Neuschae- 
fer et al. 1997, and references therein). Le Fevre et al. (1999) 
found a similar, better constrained result with p = 2.7 ± 0.6 
for < z < 1. These compare well with p = 2.2 ± 0.7 for 
Mgll absorbers with W r > 1.0 A ( |Steidel & Sargent 1992| ; the 
evolution is probably stronger for those with W r > 2.0 A). 

As such, galaxy pair evolution remains a plausible scenario 
as a contribution to the evolution of large equivalent width, 
higher redshift Mg II absorbers and the class we loosely have 
called "Double-DLAs". We note, however, that outflows from 
very luminous, star bursting galaxies at these high redshifts 
could also account for some of these "super" systems, since the 
low ioniz ation gas can be very prominent in absorption ( Pettini 
et al. 1999). 



7. SUMMARY 

We have performed a multivariate analysis of the absorption 
strengths and kinematics of Mgll and the absorption strengths 
of Fell, Lya, and Civ for a sample of 45 Mgll absorption- 
selected systems. Descriptions of the survey and analysis of the 
data used herein have been presented in Churchill et al. (1997, 
1999a) and in Paper I. The multivariate analysis was performed 
using the STATISTICA package (www.statsoft.com). We ap- 
plied both a "Tree Clustering" and "X-means Clustering" anal- 
ysis to the 30 systems in Sample CA from the data listed in 
Table 0. 

We have also compared the low, intermediate, and high ion- 
ization absorption properties of 16 systems with the B and K 
luminosities, B — K colors, and impact parameters of their host 
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galaxies. The full complement of absorption line data was taken 
from Paper I. The galaxy properties, listed in Table ||, are taken 
from the survey of Steidel et al. (1994) and the previous studies 
of Churchill et al. (1996) and Steidel et al. (1997). We tested 
for absorption-galaxy property correlations using Kendall and 
Spearman non-parametric rank correlation indicators as imple- 
mented with the program ASURV (L aValley, Isobe, & Feigel- 
son 1992), which incorporates upper limits on the data. 
The main results are summarized as follows: 

1 . The clustering analysis revealed that there is a wide range 
of properties for Mg Il-selected absorbers and that these can be 
categorized into three main types (Figure ||). Based upon strong 
W-^Lya) and W^Fell), there is the class of DLA/Hl-Rich 
Mgll absorbers. For the remaining systems, there is a spread 
in the Civ strengths for a given W^Mgll) that gives rise to 
'Civ-weak" and "Civ-strong" systems. 

2. The Civ-weak class separates into the Single/Weak and 
the C iv-deficient classes (Figures || and |). Single/Weak sys- 
tems are characterized by a single unreso ved Mg II absorption 
line with W r (Mgll) < 0.15 A, a velocity spread of ui v < 
6 km s _1 , and a range of Civ strengths, but with W r (Clv) no 
greater than ~ 0.5 A (in our sample). The Mgll strengths and 
kinematics are the dominant properties in distinguishing them 
as a separate class. C iv-deficient systems have multiple Voigt 
profile components, N c i > 1, a range of kinematic spreads, 
15 < uj v < 45 km s~\ and W r (Clv) no stronger than 0.4 A. 
The C IV strength, and to a lesser degree the Mg II kinematics, 
are the dominant properties distinguishing this as a C iv-weak 
class. 

3. The Civ-strong class separates into the Classics and the 
Doubles (Figures ^ and |J). Classic systems, like the ClV- 
deficient systems, have multiple Voigt profile components, and 
a similar, but slightly larger, range of kinematic spread. The 
Civ strengths, on the other hand, are greater than 0.4 A. Sepa- 
rate from the Classics are the Doubles, which are characterized 
by roughly twice the number of Voigt profile components, N c i, 
twice the Civ strength, and twice the Mgll kinematic spread, 
as compared to the Classics. 

4. The Mgll, and Lya strengths are tightly correlated with 
the number of Mg II Voigt profile components, N c i (see Fig- 
ure^). The tight correlation between W^Lya) and N c i implies 
that the majority of the neutral hydrogen is arising in the phase 
giving rise to the Mgll absorption. In single cloud, N c i = 1, 
systems (the Single/Weak class), W^Lya) ranges from 0.2- 
1.2 A and W r (Clv) ranges from less than 0.15 to 0.7 A, sug- 
gesting a range of metallicities and/or ionization conditions in 
this class of object. 

5. There is a highly signific ant correlation between W r (C iv) 
and uj v (see Figure H; also see Churchill et al. 1999b). This cor- 



relation is driven by the five absorbers (two Classics and three 
Doubles) with the largest uo v . For uo v < 60 km s _1 there is 
a~ 1 A spread in W r (Clv), due to the C iv-deficient systems 
and the Single/Weak systems with larger W r (C lv)/W r (Mg II) . 
More systems with uo v > 80 km s -1 are needed in order to de- 
termine if all "kinematically active" Mg II absorbers have such 
strong Civ. 

6. Assuming the Mg II clouds are in photoionization equilib- 
rium, we showed that in many systems a substantial fraction of 
the C IV absorption is arising in a separate phase from the Mg II 
(see Figure fy. We also showed that some fraction of the sys- 
tems likely have Si IV in a separate phase from the Mgll. The 
FOS profiles of the larger Civ absorbers (Classics with large 



Civ and Doubles) are resolved due to the kinematic spread, or 
velocity structure, of the C IV phase. A quantity expressing the 
degree to which the profiles are resolved, F r , correlates at the 
99.99% confidence level with the Mgll kinematic spread. Fur- 
thermore, there is a clear impression that the asymmetries in 
the resolved C IV profiles trace the Mg II kinematics. The above 
facts lead us to suggest that, often, some of the Civ arises in 
a physically distinct phase from the Mg II gas, but is kinemati- 
cally clustered with the Mg II clouds. 

7. For our small sample of 16 galaxies, we find no signifi- 
cant trends (greater than 3 a) of Civ and Lya, or their ratios 
with each other and with Mg II, with host galaxy B and K lu- 
minosities, B — K colors, and impact parameters. DLA/Hl- 
Rich systems, which have systematically smaller impact pa- 
rameters and redder colors, significantly affect the outcome of 
correlation tests. The Civ and Lya absorption properties do 
not smoothly depend upon global galaxy properties. At 2.8 a, 
we find that W r (Silv)/W r (Clv) decreases with impact pa- 
rameter, which could be interpreted as very tentative evidence 
for aglobal galactocentric ionization/density gradient (see Fig- 
ure |llj). The available data also show marginal trends for redder 
galaxies to have lower ionization conditions (Figure [l~2|). This 
is suggestive that C IV phases may be systematically smaller in 
the redder galaxies. 

8. We found that the Mgll kinematics of a given absorber 
is related to its location in the W r (Clv)-W r (Mgll) plane. In 
a comparison of our intermediate redshift sample (0.4 < z < 
1.4) with published low redshift (0 < z < 0.6) and high red- 
shift (1.2 < z < 2.1) samples, we found evidence for redshift 
evolution in the W r (Clv)-W r (Mgll) plane. This implies that 
other "classes" of absorbers are present at high redshift, pos- 
sibly including a "Double DLA/Hl-Rich" class. There is also 
a group of "super" Mg II systems at z ~ 2 that have a large 
range of W r {C\v) values. We disussed a scenario in which the 
evolution of these strongest Mg II absorbers could be due to the 
evolution of the number of galaxy pairs and/or accreting LMC- 
like satellite galaxies. 

7.1. Further Ruminations 

In general, the observed range of Mg II kinematics and C IV 
absorption strengths could be due to a number of factors. These 
include global differences in gaseous conditions related to envi- 
ronment (i.e. group or isolated galaxies) or galaxy morphology, 
local variations dominated by passage through different parts of 
host galaxies (e.g. halos, outer disks, spiral arms, etc.), or small 
scale fluctuations in galactic interstellar and halo structures. An 
additional important issue is whether characteristic absorption 
properties might be tied to the evolutionary stage of the host 
galaxy. This would imply that a given host galaxy may not al- 
ways be of the same "absorber class" throughout its evolution. 

At higher redshifts (z > 2), galaxy-galaxy interactions were 
more common and no doubt played an important role in the 
kinematics and multiphase ionization condit ions in many Mg II 
abso rbers (e.g. |Carilli & van Gorkom 1992|; B owen & Blades 
1993; [Bowen et al. 1995 ). A local analogue is seen in absorp- 
tion associated with M61, a nearly face-on, later-type galaxy 
with enhanced star formation. The Mgll velocity spread is 
~ 300 km s _1 at a n impact parameter 21 kpc and there are 
two nearby galaxies (Bowen et al. 1996). 

Since it is the kinematic complexity and spread and not the 
total amount of the Mg II that appears to be the important prop- 
erty for predicting the C IV strength in z ~ 1 galaxies, it is not 
unreasonable to conjecture that sources of mechanical energy 
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(i.e. active star formation) m ay be of central importance (see 
Guillemin & Bergeron 1997). Extended multiphase halos are 
seen preferentially around late-t ype spiral gala xies with signs 
of elevated star formatiopn (see Dahlem 1998 and references 
therein). In fact, a basic scenario which might explain both the 
taxonomy of Mg II absorbers (especially the Double and C IV- 
deficient systems) and the ClV-Mgll kinematics connection, 
is one in which star formation processes accelerate hot, ion- 
ized gas outward, which subseque ntly cools and fragments i nto 
parcels of low i onization gas (e.g. |Houck & Bregman 1990|; Li 
& Ikeuchi 1992; |Avillez 1999[ ). These parcels would be ob- 
served as outlying kinematic components in the Mgll profiles. 

These points are somewhat suggestive of a causal connection 
between evolution in the cosmic star formation rate and absorb- 
ing gas cross sections, kinematics, and ionization conditions. 
If true, we could hypothesize that scenarios of a cosmic mean 
history of galaxy evolution inferred from the global record of 
star formation history should be fully consistent with one based 
upon an "absorption line perspective". 

The global star formation rate appears to drop belo w z ~ 1.0 
( Lilly et al. 1996 ; Madau, Lucia, & Dickinson 1998 ; however, 
see |Cowie, Songaila, & Barger 1999[ ), implying" a reduction 
in the (infalling) gas reservoir for galaxies. It could be that 
a majority of galaxies transform into self-regulatory systems 
by z ~ 1 and begin to evolve in a more isolated fashion and 
in a direction dependent upon their ability to continue form- 
ing stars. Galaxies that are rich in gas and capable of forming 
large numbers of molecular clouds in their interstellar media 
(i.e. late-type galaxies) would continue to form stars and ex- 
hibit evolution, whereas those less capable would exhibit no 



discernable evolution. Such a scenario is consistent with the 
differential luminosity evolution reported by Lilly et al. (1995) 
and is not inconsistent with the tentative trends we see between 
galaxy color a nd ionization conditions (also see G uillemin & 
Bergeron 1997). 

It is unfortunate that there currently are no Mg II absorption- 
line samples for the highest redshifts (z > 2.2), so an absorp- 
tion line perspective cannot yet be fully appreciated [however, 
see the theoretical predictions of Rauch, Haehnelt, & Steinmetz 
(1997)]. 
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a Not "flag" transition; see system notes (§3) 

b Measurement of lu v was censored (see text). Uncensored values are 46.8 km s _1 (0117 + 213, z = 0.7291), 45.1 km s _1 (0117 + 213, z = 1.3430), 
c Sample A: systems with ^(Mgll) detection limits greater then 0.3 A. Sample CA: systems used for multivariate Cluster Analysis 
Sample IC: systems used for Ionization Condistions analysis. 
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TABLE 2 
Galaxy Properties 
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This paper 
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0.7291 


-22.3 


-26.3 


4.0 


36.0 


Churchill etal. (1996) 


0454 + 039 


0.8596 


-19.5 


> -23.4 


< 3.9 


10.6 


Churchill etal. (1996) 




1.1532 6 


-21.0 


-23.8 


2.8 


40.1 


This paper 
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This paper 
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Steidel etal. (1997) 
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Churchill et al. (1996), Steidel et al. (1997) 
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15.0 


Churchill et al. (1996), Steidel et al. (1997) 


2128- 123 
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3.3 


31.6 


This paper 
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This paper 



a Recent analysis of a WFPC2/HST image revealed a previously unseen galaxy within 1" of the QSO. This 
galaxy identification should be viewed with caution. 

b Mb extrapolated from 5R-band image using an Im spectroscopic template consistent with the R — K color. 
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Mgn and Civ Equivalent Widths" from the Literature 



QSO 


%abs 


W(Mgll) [A] 


W(Ctv) [A] 


References 


Mrk 205 


0.004 


0.29 


0.19 


BB93 


1219 + 047 


0.005 


1.96 


1.28 


BBP96 


0955 + 326 


0.005 


2.56 


0.85 


B99,J98 


1543 + 489 


0.075 


0.61 


< 0.26 


BBP95 


1137 + 660 


0.1164 


0.50 


< 0.32 


B94 


1704 + 608 


0.2216 


0.45 


0.43 


B94 


1317+277 


0.2891 


0.33 


0.45 


B94 


1634 + 706 


0.6694 


0.29 


< 0.07 


B94 


0454 + 034 


1.1535 


0.57 


0.95 


SS92 


1101 - 264 


1.1875 


0.41 


0.97 


BB85 


1101 - 264 


1.2030 


0.59 


0.57 


BB85 


1247+ 265 


1.2233 


0.48 


0.67 


SS92 


1602 - 002 


1.3245 


0.64 


< 0.43 


LTW87 


0226 - 035 


1.3277 


0.83 


0.43 


SS92 


1421 + 122 


1.3605 


0.43 


0.34 


LTW87 


0237 - 232 


1.3647 


2.05 


1.52 


SS92 


0957+ 561 


1.3901 


2.10 


0.27 


BB84 


0058 + 016 


1.4641 


0.45 


0.40 


SS92 


0002 - 422 


1.5413 


0.48 


0.71 


LTW87 


0424 - 131 


1.5615 


0.38 


1.05 


SS92 


1329 + 412 


1.6012 


0.70 


1.38 


SS92 


1017 + 276 


1.6081 


0.11 


0.37 


SS92 


0237 - 232 


1.6363 


0.78 


0.08 


SS92 


0421 + 016 


1.6375 


0.35 


1.13 


SS92 


1246 - 054 


1.6466 


0.52 


0.79 


SS92 


0237 - 232 


1.6583 


0.76 


1.14 


SS92 


0237 - 232 


1.6732 


1.31 


1.63 


SS92 


1311 - 270 


1.6860 


0.83 


< 0.38 


LTW87 


0958 + 551 


1.7320 


0.37 


1.01 


SS92 


0854+ 191 6 


1.7342 


4.14 


1.14 


SS92 


1017 + 276 


1.7953 


2.04 


1.51 


SS92 


1331 + 170 


1.7860 


1.05 


0.17 


SS92 


1331 + 170 


1.7766 


1.33 


1.63 


SS92 


ai AA i 1 OA 

U1UU + lol) 


1.7971 


1.11 


A AO 

0.9o 




1 TA A i 71 H 

1 f(J4 + 71U 


1 1 AT 

1.810( 


A C O 

0.58 


A C A 

0.54 


C C AO 


1 OOA i A 1 O 


1.8,555 


A A A 

0.49 


A 1 O 

O.lo 


C CAT 


1101 — 264 


1.8386 


1.00 


0.31 


LTW87 


1228 + 077 


1.8966 


1.72 


0.90 


LTW87 


1435 + 635 


1.9235 


1.05 


0.46 


SS92 


0551 - 366 b 


1.9613 


5.80 


2.13 


LTW87 


0119-044 


1.9638 


0.52 


1.68 


SS92 


0013 - 003 6 


1.9713 


7.37 


2.19 


SS92 


2126-158 


2.0219 


0.67 


0.99 


LTW87 


0013 - 003 6 


2.0290 


5.75 


0.63 


SS92 


0348 + 061 


2.0237 


0.59 


0.82 


SS92 


0424 - 131 


2.0344 


1.10 


< 0.12 


SS92 


0450- 131 


2.0667 


2.20 


1.07 


SS92 



a Equivalent widths are rest frame. 
6 Does not appear on Figure [IJ. 

References:— (BB84) Bergeron & Boisse 1984; (BB85) Boisse & Bergeron 1985; 
(LTW87) Lanzetta et al. 1987; (BB93) Bowen & Blades 1993; (B94) Bergeron 
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Fig. 1.— Three dimensional scatter plots of (a) W r (Mgll) vs. W r (Lya) vs. W r (Clv), (b) W r (Mgll) vs. W r (Lya) vs. W r (Fell), 
(c) W r (Mgll) vs. cj v vs. W r (Clv). Upper limits on Fell and Civ, of which a few are present near the origin, are plotted but not 
delineated. Vertical lines intersect the x-y values (right handed coordinate system) in the x-y plane. The z value of a datum can be 
obtained by projecting a line from the point to the background mesh in a direction parallel to the constant z lines in the mesh. 
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Fig. 2. — Tree cluster diagram showing linkage distance in a subsample of 30 absorbers, with each labeled across the bottom. The 
discriminating absorption properties are W r (Mgll), W r (Lya), W r (Clv), W r (Fell), and oj v , using a iV(0, 1) standardization (see 
text). A Ward's amalgamation algorithm was used with Euclidean distance. The isolated branches with linkages greater than ~ 4 are 
given classifications based upon the cluster means (see Figure 0). 
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Fig. 3. — .ftT-means cluster diagram for K — 5 showing the mean values W r (Mgll), Wr(bya), W r (Clv), W r (Fell), and u> v , placed 
on a N(0, 1) standardization (see text), for each of the five clusters shown in Figure With this standardization, zero is the mean 
value and ±1 is the standard deviation for a given absorption property. The cluster classifications are given on the left hand side. 
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Fig. 4. — The rest-frame equivalent widths of the elemental species listed in Table |]and in Tables 3 and 4 from Paper I vs. that 
of Mgll A2796. Also included is Call. The panels are in order of increasing ionization potential from the upper left to the lower 
right. Data point types denote the various absorber classes, and are the same as in Figure 0. Downward pointing arrows denote upper 
limits. 
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Fig. 5. — The rest-frame equivalent widths of Mgll, Fell, Ly a, and C IV vs. the Mgll kinematic spread, u> v . Data point types denote 
the various absorber classes, and are the same as in Figure pi 
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Fig. 6. — The rest-frame equivalent widths of Mgll, Fell, Lya, and C IV vs. the number of Voigt profile components, or clouds, N c i. 
Only Sample IC is presented for reasons given in the text. Data point types denote the various absorber classes, and are the same as 
in Figure H. 
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Fig. 7. — The maximum W r (Clv) (left panel) and W r (Silv), under the assumption of extremely high photoionization equilibrium, 
that can arise in single phase Mgll clouds vs. the measured rest-frame equivalent widths. The predicted points are based upon 
CLOUDY models of the systems (Sample IC), incorporating the number of clouds and their column densities, b parameters, and 
kinematics (see text for details). If a system lies below the line, then the measured Civ or Si IV cannot arise solely in the Mgll 
clouds, but must arise in a distinct highly ionized phase that does not give rise to detectable amounts of Mg II. The typical error bar 
of the measure data are given in the upper left corners of both panels. Data point types denote the various absorber classes. 
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Fig. 8. — The Civ kinematics. Shown are the Mgll A2796 transitions and Civ AA1548, 1550 doublets of the 17 systems for which 
both members of the Civ doublet are detected and unblended with other transitions in the FOS spectra. The panels with Mgll show 
a velocity window of 600 km s _1 , centered on the system zero point. Ticks above the continuum give the Voigt profile component 
velocities and the smooth curves through the data are the Voigt profile fit results. The panels with the Civ doublets show a velocity 
window of 2000 km s _1 . The Mgll Voigt profile velocities are projected over the continuum of each doublet. The smooth curves 
through the data are single Gaussian fits with the widths held at the instrumental resolution, but with the centroids and depths allowed 
to vary. 
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Fig. 9. — The Fell and Civ kinematics of a DLA at z = 1.7763 toward Q 1331 + 170; the Mgll profiles were not covered due to the 
setting of the HIRES echelle. The upper panels show HIRES/Keck data at R = 22, 500 over a velocity window of 500 km s _1 . Ticks 
above the Civ roughly give the kinematic composition of the high ionization gas. The lower panel shows simulated FOS spectra of 
the Civ doublet over a velocity window of 2000 km s _1 . As with Figure || single Gaussian fits with widths held at the instrumental 
resolution are superimposed on the data. 
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Fig. 10.— The Mgll kinematic spread, andMgll, Fell, Lya, and C IV absorption strengths vs. host galaxy impact parameter, absolute 
K magnitude, and rest-frame B — K color. 
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Fig. 11.— Selected equivalent width ratios vs. host galaxy impact parameters. 
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Fig. 12.— Selected equivalent width ratios vs. host galaxy B — K colors. 
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Fig. 14. — The W r (Clv)— Wr(Mgll) plane. Data from this study have point types the same as in Figure |[ Also plotted are: higher 
redshift points (1.2 < z < 2.2) from Bergeron & Boisse 1984, Boisse & Bergeron 1985, Lanzetta et al. (1987) Steidel & Sargent 
(1992) [six-pointed stars], lower redshift points (0.1 < z < 0.6) from Bergeron et al. (1994) [downward pointing, open triangles], 
and z ~ data from Bowen et al. (1996, 1996), and Bowen (private communication, 1999) [open diamonds]. 
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Fig. 15. — The Mgll A2796 profiles for four higher redshift systems and for the z ~ SN 1993J, a line of sight from the galaxy M81. 
These systems all have equivalent widths greater than 2 A. They occupy the extreme right portion of the C iV-Mg II plane (Figure |l4j) 
and are members of the population of the most rapidly evolving Mgll absorbers. 



